Astronomy & Astrophysics manuscript no. AstroPHpaper 


©ESO 2012 


September 11, 2012 





f 



The reawakening of the sleeping X-ray pulsar XTE J1 946+274 

Sebastian Miiller 1 , Matthias Kiihnel 1 , Isabel Caballero 2 , Katja Pottschmidt 3 , Felix Fiirst 1 , Ingo Kreykenbohm 1 , 
Macarena Sagredo 14 , Maria Obst 1 , Join Wilms 1 , Carlo Ferrigno 5 , Richard E. Rothschild 6 , and Riidiger Staubert 7 

1 Dr. Karl Remeis-Observatory & ECAP, Universitat Erlangen Niirnberg, Sternwartstr. 7, 96049 Bamberg, Germany 

2 CEA Saclay, DSM/IRFU/SAp-UMR AIM (7158) CNRS/CEA/Universite Paris 7, Diderot, 91191 Gif sur Yvette, France 

3 CRESST and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA and Center for Space Science and Technology, 
UMBC, Baltimore, MD 21250, USA 

4 Department of Physics, Florida International University, Miami, FL 33199, USA 

5 ISDC Data Center for Astrophysics, University of Geneva, 16 Chemin d'Ecogia, 1290 Versoix, Switzerland 

6 Center for Astronomy and Space Sciences, University of California, San Diego, La Jolla, CA 92093, USA 

7 Institut fur Astronomie und Astrophysik, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 



Received / Accepted 



ABSTRACT 



f 



We report on a series of outbursts of the high mass X-ray binary XTE J1946+274 in 2010/201 1 as observed with INTEGRAL, RXTE, 
and Swift. We discuss possible mechanisms resulting in the extraordinary outburst behavior of this source. The X-ray spectra can be 
described by standard phenomenological models, enhanced by an absorption feature of unknown origin at about lOkeV and a narrow 
iron Ka fluorescence line at 6.4 keV, which are variable in flux and pulse phase. We find possible evidence for the presence of a 
cyclotron resonance scattering feature at about 25 keV at the 93% level. The presence of a strong cyclotron line at 35 keV seen in data 
from the source's 1998 outburst and confirmed by a reanalysis of these data can be excluded. This result indicates that the cyclotron 
line feature in XTE J1946+274 is variable between individual outbursts. 

Key words. X-rays: binaries - pulsars: individual XTE J1946+274— Accretion, accretion disks 



1. Introduction 

Due to the ~10 12 G strong magnetic field at the magnetic poles 
of many accreting neutron stars in high-mass X-ray binaries, cy- 
clotron resonance scattering features (CRSFs or cyclotron lines) 
are observable in the X-ray spectra of these sources. These lines 
originate from photons generated in the accretion column of a 
magnetized neutron star interacting with electrons in the column, 
since their motion perpendicular to the B-field is quantized into 
Landau-levels with energy differences 



AE * 12keV 



(io 12 g) 



(1) 



To date CRSFs have been reported for about 20 X-ray pulsars 
(ICaballero & Wilmsll201 lh . 

The 15.8 s pulsar XTE J 1946+274 was first detected in 1998 



(ISmith & Takeshimal[l998l: IWilson etal 



19981). It is a transien t 



X-ray source with a Be-type companion ( Verrecchia et alj2002h . 
This kind of main-sequence B stars shows Balmer emission lines 
caus ed by a quasi-Keplerian equatorial disk near the Be star (see, 
e.g. JHanuschlll99ftlSlettebaklll988h . Interaction between this 
disk and the orbiting neutron star can lead to violent X-ray out- 
bursts resulting in the appearance of a bright X-ray sources in 
the sky. 

The initial outburst of XTE J 1 946+274 in 1998 lasted about 
three months. iHeindl et al.l d200ll) reported on the discovery of 
a CRSF at an energy near 35keV during this outburst, which 
was followed b y sev eral fainter outbursts separated by ~80d 
(ICampana et alJll999h . This separation was later established by 
IWilson et al.l (120031) as half of the ~170d orbital period. The oc- 
currence of two outbursts per orbit could be related to the mis- 
alignment of the Be star's angular momentum and the orbital 



plane of the neutron star. While we look n early onto the pole of 
the Be star, the orbital inclination is >46° (IWilson et alj|2003l) . 

After a phase of nearly periodic flaring between 1998 and 
2001, XTE J1946+274 went into quiescence unti l 2010 June 
when a new sequence of outbursts started (see, e.g jMuller et al.l 
120101) . The principal outburst behavior in 2010 was the same 
as that observed in the 1990s, with two outbursts per orbit. As 
shown in Fig. Q] however, the five outbursts seen during 2010 
are neither clearly connected to the times of periast ron and apas- 
tron p assages (based on the orbital ephemeris from lWilson et al.l 
120031) . nor to any other specific orbital phase. This new outburst 
episode ended in 2011 June and the source was back in quies- 
cence at the time of writing. 

In this paper, we report on INTEGRAL, Swift, and RXTE 
observations of the first and the third outburst of this series in 
2010 June/July and 2010 November/December, respectively. In 
Sect. [2] we give a summary of the observations. In Sect.[3]we de- 
scribe the continuum model and study the time and pulse phase 
resolved behavior of the spectra. Summary and conclusions are 
given in Sect. [4] 

2. Observations and data reduction 

The 2010/2011 outburst series started on 2010 June 4, 
when SwiftfBAT detected an increase of the X-ray flux of 
XTE J194 6+274, rising up to 40mCrab (15-50keV) within 
three days dKrimm et al. I l20l0h . This first outburst lasted about 
60 d and reached a flux of ~140mCrab (see Fig. [TJ. The subse- 
quent four outbursts lasted between 30 and 50 d, each, reaching 
almost the same maximum flux level between 40 and 60 mCrab. 
The separation of the outbursts is between 60 and 90 d. A pecu- 
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Table 1. Summary of all observations used. 
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Fig.l. 15-50keV SwiftfB AT light curve of the 2010/2011 out- 
burst series. The times of periastron and apastron passages are 
marked by P and A, respectively. These epochs and the corre- 
spon ding uncertainties we re calculated using the orbital solution 
from IWilson et alJ (120031) . The insets provide a closer view on 
the first (Jun/Jul) and the third (Nov/Dec) outburst. The blue, red, 
and green arrows in these insets indicate the observations times 
for PCA, ISGRI, and XRT, respectively. Epochs over which data 
were summed for the time resolved spectral analysis are indi- 
cated with I to VI. 



liar behavior was observed during the third outburst where the 
primary maximum was followed by another brightening, reach- 
ing again luminosities up to ~50mCrab (Fig.[TJ inset). 

We present data from the Proportional Counter Array 
(PCA, Uahoda et al l 120061) on board the Rossi X-ray Timing 
Explorer (RXTE. iBradt et alJ 1 19931), the INTEG RAL Soft 
Gamma-Ray Imager (ISGRI, iLebrun et alJ 120031) on board 
the INTErn ational Gamma-Ray Astrophysics Laboratory 
(INTE GRAL. IWinkler et al l 120031) . and the X-ray Telescope 
(XRT, iBurrows et al.l I2Q05D on board the Swift Gamma-Ray 
Burst Explorer ( Gehrelsl 20041) . Data were reduced with the 
standard analysis pipelines, based on heasoft (v. 6.10 and 6.1 1) 
and INTEGRAL OSA v. 9.0. XTE J1946+274 was monitored 
by RXTE regularly during the first and the third outburst. The 
source was also sporadically in the field of view of INTEGRAL 
during the first and the third outburst. Swift pointings were 
available only during the main and secondary peak of the third 
outburst. Table Q] contains a log of the observations with these 
satellites, which are also indicated in Fig.Q] 

PCA consisted of five proportional counter units (PCUs) 
with a field of view of ~1°, sensitive between 2 and 90keV. 
Since PCU2 is known to be the best calibrated one dJahoda et alJ 
120061) . only data from the top layer of this PCU are used. We ob- 
tained 2-60 keV light curves with 0.125 s resolution, spectra in 
the standard2£ mode, and pulse phase resolved spectra using 
GoodXenon data. The light curves were corrected to the barycen- 
ter of the solar system using faxbarjfl. The PCA background 
model SkyVLE was used for PCA background sub traction. Due 
to lar ge uncertainties in the orbital parameters (IWilson et alJ 
120031) no correction could be performed for the neutron star's 
orbital motion. The lack of sufficient statistics prevented us from 
improving the existing orbital solution. Data from the High 
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Notes. {a) For PCA, the first column contains the number of the Obs-ID 
after 95032-12- and for XRT after 0031888. For ISGRI the revolution 
number is listed. {b) Total background corrected counts. Negative count 
rates are due to uncertainties in the background determination. (c) Epoch 
for data grouping, see text for details. ObsIDs without assignment in 
this column are different enough to other spectra during the same epoch 
that they cannot be combined with the other data or they were recorded 
without simultaneous low-energy measurements. We therefore exclude 
these data from our analysis. 



http://heasarc.gsfc.nasa.gov/lheasoft/ftools/fhelp/faxbary.txt 



Energ y X-ray Timing Experiment (HEXTE, iRothschild et al.l 
1 19981) on board RXTE were excluded from our analysis since 
both HEXTE clusters were not rocking at the time of the obser- 
vations. The resulting uncertainties in the background determi- 
nation are too large for the purposes of this paper. 

The CdTe detector of INTEG RALflSGRl covers th e energy 
range from ~18keV to 1 MeV dLebrun et al.1 120031) . Thanks 
to the large field of view of INTEGRAL, XTE J1946+274 
was detected several times during its recent outbursts between 
2010 June and 2011 April in observations pointed at Cyg X-l. 
We extracted ISGRI pulse phase averaged spectra for all obser- 
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Fig. 2. Panels a and b show two example spectra of 
XTE J1946+274. Epoch I is the spectrum with the highest coun- 
trate while epoch IV is the only data set where all instruments are 
available. Histograms show the best fit continuum. Panels c-k 
show the behavior of the residuals when adding the additional 
components of the spectral model one by one. Numbers at the 
bottom right of each panel indicate the best-fit statistics, ^- 2 /dof. 
c and d: best fit using the Fermi-Dirac cutoff only, e and f: resid- 
uals after adding the lOkeV feature to the continuum model, g 
and h: residuals after adding the Galactic ridge emission to PCA 
(for epoch I the flux of this model component was fixed to the 
result obtained from epochs III and IV). i and j: residuals after 
adding the Fe Ka line to the model. Since in the model fit of 
epoch IV the depth of the CRSF results to zero, panel j displays 
the final result for epoch IV. In panel k, the CRSF has been added 
to the data for epoch I. See also text for a discussion of the sta- 
tistical significances of the Galactic ridge emission, the lOkeV 
feature, and the CRSF. 



vations of XTE J 1 946 +274 using the standard spectral extraction 
method of OSA 9 described in the INTEGRAL documentation^. 
We selected those observations of XTE J1946+274 for which it 
was less than 12° off-axis. 

To cover the soft X-ray band we use data from Swift/XKT, a 
600 x 600 pixel CCD covering a field of view of 23 '.6 x 23. '6 in 
the energy range 0.2-10keV. We extracted the data in Windowed 
Timing mode. For the source region we have choosen a circle 
around the source with a radius of ~V . The background was 
extracted using circles of the same radius at two off-source posi- 
tions about 2'5 away from the source. 

In the spectral analysis, for XRT we use data between 
1.5 keV and 7.0 keV. We binned the XRT and PCA data to a 
signal-to-noise ratio of ten. Due to calibration problems of PCA 
at the Xe L-edge, we exclude PCA data below 4.5 keV from 
our analysis. Discarding bins with S /N < 10 results in an up- 
per energy limit of typically 40 keV for the PCA. ISGRI covers 
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Fig. 3. Results of time resolved spectroscopy. The units are 
as follows: Fx in 10~ 2 photons s cirT 2 in the 7-15 keV en- 
ergy band, Nh in 10 22 cirT 2 , Ef \a and E cut in keV, and Ap e in 
10~ 4 photons s~' cm~ 2 . 



the high energy band between 20keV and 100 keV. We added a 
systematic error in quadrature to the PCA and the ISGRI spec- 
tra using canonica l values of 0.5% and 2.0%, respectively (see, 
iJahoda et al]|2006l and IBIS Analysis User Manual). 



3. Spectral analysis 

We performed all fits usi ng the Interactiv e Spe ctral 
Interpretation System (ISIS, iHouck & Denicolal l2000l) . In 
order to improve the signal to noise ratio of individual spectra 
we averaged the data over six data blocks in time taking into 
account the flux level and instrument availability. These epochs 
I- VI are defined in Table [2] and shown in Fig. Q] The spectra 
in epochs I and II cover the first outburst at the maximum and 
fading phase, respectively. Epochs III-V follow the maximum 
and the fading phase of the third outburst. The flux in the 
maximum level of this outburst (epoch III) is comparable to 
the fading phase of the first outburst (epoch II). The last set of 
observations is summarized in epoch VI. These XRT data cover 
the fading phase of the flare right after the third outburst. Before 
defining these epochs we confirmed that the spectral variability 
during these epochs is negligible. 



3. 1 . Spectral model 

As shown, e.g., bv lBecker & Woil J2007T) . the X-ray spectra of 
accretion-powered X-ray pulsars can be roughly described by a 
powerlaw with a high energy cutoff. In practical data modelling, 
this continuum has been approximated by several di fferent con- 
tinuum models (see, e.g.. iRrevkenbohm et alj|2002l for a sum- 
mary). Here w e describe the data using the so called Fermi-Dirac 
cutoff (FDCO. lTanakalll986l) . given by 



FDCO(£) oc E x 



1 + exp 



E E cl 

£fold 



(2) 



http://www.isdc. unige.ch/integral/analysis#Documentation 



which has been successf ully applied to other accre ting X-ray 
pulsars such as Vela X-l dKrevkenbohm et ai]|2008l) . This con- 
tinuum is modified by the CRSF, modeled as a line with a 
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Table 2. Results of the time resolved spectral analysis. 



e" 


start 6 


stop 6 


cts c 


Fi d 


F 2 d 




r 


-g — j 


~E~t 


A Fe « 




X 1 ,/dof 

A mil' 


I 


367 


378 


17.20 


6.91(2) 


8.40(2) 


2 ii-, 

—2 1 


75 + ™ 

" -0.04 


7 25 + »'^ 

-0.30 


16.71^ 


10 ± 4 


55 ± 20 


1.30/66 


11 


380 


394 


16.46 


4 71(3) 

-r. / j. \~> / 


5.72(2) 


2 6 + li 

—2 5 


82+fH 4 , 

—0 32 


6 1 + \l 

-1.0 


19 4 +2: i 

-9.7 


5 9 + li 

^ -2.5 


46 ± 20 


1.03/45 


III 


524 


529 


2.93 


3 03(31 


3.70(2) 


| 77+O.25 
l. / '_ 2 9 


74 +o:l2 

'^-0.17 


8 +2 4 


16 +4 


1.5 + 1.0 


18 + 14 


1.17/202 


IV 


530 


534 


1.68 


2.54(3) 


3.17(2) 


L84+ 24 


°- 74+ o?o 


8- 1+ 06 


14 + 4 


°- 4+ 04 

— U.4 


6+' 5 


1.03/176 


V 


535 


538 


1.12 


1.88(2) 


2.46(2) 


5 1+ 2 - 5 ' 




6-01?! 


19.0!^ 


< 1.5 


< 30 


0.78/32 


VI 


552 


559 






1.69(10) 


7 14+025 


0.84!«»f 


7.25 


16.7 


< 1.0 


< 27 


1.21/135 






c pca' 


t'lSGRl' 


c xrt' 


E G > 


t g j 


<T G } 


W " 


^CRSF^ 


T CRSF ; 


O" CRSF^ 


WcRSF ' 


I 


0.943(20) 




0.97(3) 




9 85 +u - 2U 
9 9+0-4 


6.9+ 17 


9 9+0.8 

-9-5 


-380 + ™ 


9<: q+0.« 


9+io 


65 +1 46 


1H "-120 


II 


1.00(7) 








3 9 +4: ^ 


1.8!™ 


_ 1 70+130 


25.3 


1+1 


0.65 


-60 +6 ° 


III 


0.97(6) 






0.88(2) 


9.8^| 
9 8+^ 


t q+il 


2 


-iso:f 


25.3 


<23 


0.65 


> -330 


IV 


0.97(4) 




0.88(7) 


0.92(2) 


lr® 

J -'-2.3 


2 


-170!™ 


25.3 


<40 


0.65 


> -500 



VI 1 — - 

VI — — — 1 
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Gaussian optical depth profile 
exp(-T(£)), with t(E) = t C rsf X exp 



1 / E - £crsf \ 

2 \ CT CRSF / 



(3) 



The equivalent width of this feature will be denoted with Wcrsf 
in the following. The spectra of some X-ray pulsars, including 
XTE J 1946+274, also contai n an absorptio n or emission like 
feature in the range 8-12keV dCoburnll200lh . The origin of this 
so-called 10 keV feature is still unclear, however, since it appears 
always at about the same energy, it is probably not related to the 
magnetic field strength of the neutron star. We modelled this fea- 
ture as a broad Gaussian absorption feature, as we also did for 
the CRSF (Eq. [3]l with the centroid energy, Eg, the width, <r G , 
optical depth, tg, and equivalent width, Wg- Finally, interstel- 
lar absorption was m odelled with an upda ted version ofTBabsQ 
using abundances by I Wilms et al ] (120001) and cross sections by 
IVerner & YakovleviJ^gsl) - 

In modelling the data, we have to take into account that the 
cross-normalization of the different instruments used is not per- 
fectly well known and the source might also slightly change in 
flux between the different observations. These effects were taken 
into account by cross calibration constants cxrt and cisgri, us- 
ing the PCA as the reference instrument. Furthermore, the PCA 
background was allowed to vary slightly in count rate. To ac- 
count for this imperfections in the modelling of the background, 
we introduced the constant Ch. 

Finally, the data modelling is affected by the proximity of 
the source to the pl ane of the Galaxy. Galactic ridge emis- 
sion (GRE, see, e.g.. IWorrall etaT1ll982t IWarwick et al.lll985l) 
manifests itself through the presence of an emission feature at 
~6-7 keV in the PCA spectrum caused by unresolved Fe Ka flu- 
orescence lines. These lines are not present in the XRT data in 
epochs III and IV and therefore must be due to diffuse emis- 
sion that is picked up by the PCA only, due to its larger field 
of view. In our modelling of the PCA data we therefore intro- 
du ced a model for the Galactic ridge emission based on that 
of lEbisawa et al ] (120071) . who described the ridge emission as 
the sum of two bremsstrahlung components and an iron line 



see http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ 



complex modelled by three Gaussian lines at 6.4 keV, 6.67 keV, 
and 7.0keV, with equivalent width ratios of 85:458:129, re- 
spectively. Since we used the PCA data down to 4.5 keV only, 
we did not account for the two soft bremsstrahlung compo- 
nents and modelled the Galactic ridge emission as the sum 
of three narrow Gaussians, with fixed energies and equiva- 
lent width ratios according to lEbisawa et al.l d2007t) . This com- 
ponent was applied only to the PCA data and absorbed by 
the interstellar medium using the Galactic Nh as determined 
from the Leiden/ Arge ntine/Bonn (LAB) Su rvey of Galactic H i, 
N H = 9.4 x 10 21 cm- 2 dKalberla et alJl2005l) . We determined the 
flux of the Galactic ridge emission from simultaneous fits to 
the XRT and PCA data from epochs III and IV. In epoch IV, 
residuals around 6-7 keV are only visible in PCA, while in 
XRT neither residuals from this emission, nor from a source 
intrinsic iron Ka line are detected. Thus, the residuals in the 
epoch IV PCA spectrum must be caused by Galactic ridge 
emission. In epoch III, on the other hand, we find weak ev- 
idence for the presence of a source intrinsic iron Ka line in 
the XRT spectrum. The unabsorbed flux of the 6.4 keV iron 
line of the Galactic ridge emission in both spectra is consis- 
tent with each other ((7.2 + 2.0) x 10~ 5 photons s~' cuT 2 , and 
(6.6 + 2.0) x 10~ 5 photons s~' cm 2 for epoch III and IV, respec- 
tively). While these values are sl ightly higher than th e 1.22 X 
10~ 5 photons s _1 cm -2 reported bv lEbisawa et alJ d2007l) . the dif- 
ference is still within the typical variation of the Galact ic ridge 
emission over the Galactic plane dYamauchi et al.ll2009l) . Since 
the flux of the Galactic ridge emission is constant over time, 
we added this model component with parameters fixed to the 
mean value as obtained from epochs III and IV to the PCA spec- 
tra of all epochs. We find that for the late part of the outburst, 
the ridge contributes 1.5% of the 3-10keV flux and 8% in the 
Fe band (6-7keV). To estimate the significance of the Galactic 
ridge emission, we performed Monte Carlo simulations of the 
best fit model without this feature to create a set of 1000 syn- 
thetic spectra. We then performed the fit allowing all model pa- 
rameters, including the Galactic ridge emission, to vary. For both 
epochs III and IV none of these simulations led to a fake spec- 
trum for which the improvement in x 1 was as large as in the 
real data, i.e., the probability that the Galactic ridge emission is 
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Fig. 4. a: PCA spectra of time bins I-V (shifted in y-direction 
for better visibility), b: Ratio of the background corrected time 
resolved spectra. 

real is greater than 99.9 % (>3.3cr). The remaining residuals at 
6.4 keV can be explained by a narrow source intrinsic iron Ka 
flourescence line. We model this feature by a thin Gaussian emis- 
sion line with fixed centroid energy Ep e = 6.4 keV and width 
cr Fe = l(T 4 keV. The flux, Ap e , (and thus equivalent width, Wp e ) 
were allowed to vary. 

In summary, the model used can be written as 

M = TBabs x (FDCO + Fe 6 . 4keV ) x Giokev X Gcrsf + GRE. (4) 



3.2. Time resolved spectroscopy 

In this Section we describe the time resolved behavior of the 
spectral parameters. For each of the six epochs we fitted the re- 
spective spectra separately, including the Galactic ridge emis- 
sion as a constant component as discussed above. Example spec- 
tra of two epochs together with the best model fit are shown in 
Fig. |2] The free fit parameters are summarized in Table [2] and 
displayed in Fig. [3] We calculated fluxes in the energy band 
10-20keV for the epochs I-V. In order to be able to compare 
the source fluxes of all epochs, but avoid excessive extrapola- 
tion of the models, we also derived the fluxes in the energy band 
7-15 keV. We excluded the contribution from the Galactic ridge 
emission from this flux. 

The resulting xLa f° r an nts does not exceed 1.3 (see 
Table |2J. For epoch I, the value of x^ ei is rather high. Here, 
the greatest contribution to originates from residuals of the 
ISGRI data which are caused by calibration problems, so we ac- 
cept this fit. Note that not all model components are necessary to 
describe the data in the spectra with low statistics, e.g., epoch V 
provides an statistically too low^ ed of 0.78. The reason for in- 
cluding these components in these fits as well, even though the 
components overdetermine the fit model, is that this way it can 
be shown that these spectra are at least consistent with the full 




0.2 0.4 0.6 0.8 
pulse phase (p 

Fig. 5. PCA-Pulse profile of epoch I (in the full PCA energy 
band, i.e., 2-60 keV), shown twice for clarity. The count rate in 
each bin is normalized to its variance relative to the mean pulse 
profile count rate. The phase bins are marked as follows. MP: 
main peak, SP: secondary peak, off: off state. 



model. In addition, we note that an overestimation of the sys- 
tematic error would also yield a Xied that is too low. In epoch VI 
only XRT data are available and the continuum parameters are 
badly constrained from these data alone. We therefore fixed Ef n 
and E cat to the value obtained from epoch I because a change of 
these parameters affects mainly energies not covered by XRT. 
However, possible influences of these fixed parameters to the 
free fit parameters A^h and F might affect the results. 

The behavior of the photon index, F, the cutoff energy, 
E cul , and the folding energy, £f id, yields information about 
the evolution of the spectral continuum. In most cases there 
are no or only slight variations of these parameters apparent. 
However, we know from previous observations that these param- 
eters can also be significantly correlated to each other. We there- 
fore also derived a model independent illustration of the spec- 
tral changes by dividing the background subtracted PCA spectra 
from epochs II- V by the spectrum from epoch I (see Fig.[4j». The 
ratios II/I, III/I, and IV/I appear to be mainly constants, meaning 
that the variations of the continuum parameters in epochs I-IV 
are probably artificial and due to cross correlations. In contrast 
to this, ratio V/I shows a spectral softening for epoch V, caused 
by a real change of the continuum parameters. 

Except for epochs covered by XRT, the hydrogen column 
density A^h is not well determined and the best-fit parameters are 
consistent with constant A^h- In addition to the line caused by 
Galactic ridge emission, an Fe Ka fluorescence line at 6.4 keV 
is required during the first outburst and the maximum phase of 
the third outburst (epochs I— III) . The flux of this line is corre- 
lated with the X-ray flux of the source, as expected for a fluo- 
rescent line. Furthermore, for these epochs at highest luminosi- 
ties (I— III), the equivalent width Wp e stays roughly constant. In 
the fading phase of the third outburst and its subsequent flare 
(epochs IV- VI), this feature is consistent with zero, meaning that 
the additional emission at these energies can be explained by the 
Galactic ridge emission. To estimate the significance of the Fe 
Ka line, we performed similar Monte Carlo simulations as those 
done for the Galactic ridge emission. For epochs I and II, the 
probability that there is a source intrinsic Fe Ka line is greater 
than 99.9 % (>3.3<x). During epochs III and IV, where the source 
was much fainter, these simulations show that the probability for 
Fe line emission from the source is 98.5% (2.4<x) and 62.0% 
(0.9cr), i.e., here the Fe line region is dominated by Galactic 
ridge emission. 
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Fig. 6. a: PCA spectra of the three phase bins together with the 
pulse phase averaged spectrum (shifted in y-direction for better 
visibility), b: Ratio of the background corrected pulse phase re- 
solved spectra. 



We find no evidence for the centroid energy, the width, and 
the optical depth of the lOkeV feature to be variable over time. 
Due to the lack of statistics, the width of the feature, cr G , in 
epoch III and IV cannot be constrained and is therefore fixed to 
the mean value obtained from epochs I and II, i.e., cr G = 2.0 keV. 
In epoch V, also due to the lack of statistics, this feature is not re- 
quired to describe the data. We therefore omit it from our model 
in this epoch, and also from epoch VI, where our coverage exists 
only below 7 ke V. 

One of the most interesting questions is whether there is a cy - 
clotron line present as in the 1998 outburst dHeindl et al.ll200 If) . 
We find a possible CRSF at ~25keV during epochs I-IV, i.e., 
those epochs where good coverage exists above lOkeV. The 
statistics during the maximum of the first outburst (epoch I) 
are good enough to obtain CRSF parameters that are well con- 
strained. As is fairly common for CRSF fits dCoburnll200il) . the 
width and depth of the line are strongly correlated. We therefore 
set a lower limit of 0.5 keV for the width, which is comparable 
to the resolution of the PCA at these energies. Monte Carlo sim- 
ulations as described above lead to a 93% (1.81<x) probability 
that the CRSF found in epoch I is real. 

The dependence of the CRSF's parameters to the choice of 
the approach of modeling the continuum is an important point 
which also has to be discussed. The results presented here were 
obtained by first adding the lOkeV feature to the model, and 
then accounting for the CRSF. Doing this vice versa for epoch I, 
we find another minimum of xL&> leading to different cyclotron 
line parameters, i.e., £crsf = 29.6^} | keV, tcrsf = 0.34^^, 
and cr CRSF = 6.8^] ^keV The main problem for this alternative 
approach is the need of fixing parameters to certain values for 
the final fit. For example, allowing cr CRSF for the final fit to vary 
leads to an unrealistically broad CRSF which effectively models 
part of the exponential rollover and not the line. Furthermore, 
the quality of this fit is slightly worse (x 2 « 97 vs.^ 2 « 86). We 
note that there is a third solutio n in which the centr oid energy 
of the CRSF is in agreement with lHeindl et alJ (|200lh . However, 
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Table 3. Results of the pulse phase resolved spectral analysis. 



parameter 


MP 


SP 


off 




0.50-0.91 


0.13-0.34 


0.34-0.50 & 0.91-0.13 


r 

Efoid [keV] 


0.87 ± 0.05 
6.7 ± 0.4 


49+0.10 

7 Q+0.4 
'■ 7 -0.5 


rv t 1 +0 07 
v,/ A -0.09 

7.8 ± 0.5 


£ cut [kevj 

A Fe LIU J 


1 o 1+1.2 
19 ' 2 -1.5 
1 l.y ± Z.J 


11+4 

o 7+2.8 

s. i_ 29 


1 ^ n+2.4 
13.U_ 3 j 

y.j ± Z.Z 


W Fe [eV] 




46+ 15 


54+14 


E G [keV] 


q 7+0.6 


9.5 + 0.3 


9.9 ± 0.2 


T G 

W a [eV] 


034 +0 015 
-2001* 


0.080!^ 
-170 +60 
04 +() ^ 

-O.04 
U "-I40 


0.114 + 0.018 
-600 + «» 


T CRSF 


0.12 + 0.08 


0.13 + 0.09 


Wcrsf [eV] 


-180 + 90 


-210+130 


Xtjdof 


0.92/56 


1.10/55 


0.93/56 



Notes. Uncertainties and upper limits correspond to the 90% confi- 
dence level for one parameter of interest. 

(a) In units of photons s cm" 2 . The centroid energy has been fixed to 
6.4 keV, and the width to 10~ 4 keV. 



contrary to iHeindl et al.l d2001l) . this third solution has an un- 
physically broad and shallow shape and is thus not physically 
meningful. All results presented in this paper are based on the 
approach first adding the lOkeV feature, and then the CRSF. 
Finally, we note that there are two further solutions, where the 
CRSF is located at 30keV and 40keV and a line width of less 
than 1 keV. These two solutions are indeed physically mean- 
ingful, however, a Monte Carlo estimation of the significances 
of these features yields probabilities of 85% (1.44<x) and 86% 
(1.48<x), respectively, much lower than for the solution with the 
CRSF at 25 keV. Furthermore, these two solutions are supported 
only by about three data bins per instrument (the 40 keV solu- 
tion even only by ISGRI), while the 20 keV solution is based on 
a much larger number of spectral bins. 

Including a CRSF in the fainter observations in epochs II-IV 
does not significantly improve the fits. All observations are in 
principle consistent with the presence of a weak line as that seen 
in epoch I. Including for consistency such a feature in the spec- 
tral modeling, fixing the centroid energy and width and leaving 
the optical depth respectively the equivalent width of the line as 
a free parameter effectively gives an upper limit for the depth 
of the line in these observations. Not unexpectedly, the limit 
for tcrsf becomes larger for the fainter phases of the outburst 
(Table 0. We find no evidence for the equivalent width Wcrsf 
to be variable over time. 

3.3. Pulse phase resolved spectroscopy 

Spectral parameters are known to be variable as a function 
of pulse phas e, e.g., the param eters of cyclotron lines (e.g., 
1A 1 1 18-61: ISuchv et alJl201lh . To investigate such changes, 
we performed pulse phase resolved spectroscopy on epoch I 
with PCA, based on a pulse period ephemeris found by epoch 
folding. Epoch I provides the best statistics for this analy- 
sis. The ephemeris can be described with a constant spin 
up, /"puise = (-3.0 + 0.3) x 10~ 9 s s _1 with pulse period between 
15.755 s and 15.767 s. The uncertainties of the pulse periods 
were typically less than 0.003 s. Since it was not possible to cor- 
rect the light curve for the orbital motion of the system, this spin 
up trend might be caused by orbital effects and not by transfer 
of angular momentum of the accreted matter onto the neutron 
star. We extracted pulse phase resolved spectra using individual 
pulse periods following this ephemeris. The resulting pulse pro- 
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Fig. 7. Flux of the iron Ka fluorescence line against the 
7-15 keV flux. The crosses correspond to the data from the 2010 
outburst, and the diamond to the average spectrum of the 1998 
outburst, obtained by reextracting and analyzing the old PCA 
and HEXTE data from 1998. 



file (Fig. ID sho ws a double peaked s tructure, typical for such 
sources (see, e.g jBildsten et al]|1997l) . 

As in the time resolved case, we used the FDCO model to fit 
the phase resolved spectra using three phase bins as indicated in 
Fig-El The main peak (MP), the secondary peak (SP), and the off 
state (off). The pulse phases of the bins used are summarized in 
Table [3] Constant Galactic ridge emission is taken into account 
in the same way as for the pulse averaged analysis. The back- 
ground scaling factor for PCA Cb is assumed to be equal to the 
result from the pulse averaged analysis and was therefore frozen 
to this value. Initial fits of the pulse phase resolved spectra lead 
to large unconstraints in A^h and the width of the lOkeV feature 
ctq. These values were therefore frozen to the results form the 
phase averaged analysis. The signal to noise ratio in the spec- 
tra from the individual pulse phases is too low to allow us to 
study the behavior of the CRSF Similar to the flux dependent 
analysis, we freeze the CRSF centroid energy and width to the 
result from the pulse averaged analysis and allow only the op- 
tical depth respectively the equivalent width to vary (due to the 
marginal significance of the CRSF no further attempt was made 
to study a phase dependence of the cyclotron line). 

Figure [6] shows the ratios of background corrected pulse 
phase resolved spectra. The ratio MP/off appears mainly as a 
constant, while the other two ratios reveal more complex shapes. 
This constant ratio is in agreement with the spectral parameters 
of the main peak and the off state. These parameters show almost 
no or only slight differences to each other. The ratio MP/SP, a 
falling line, can be explained by a change in photon index. This 
explanation is confirmed from the spectral analysis. While E cut 
and Ef id show only small differences, the photon index F is sig- 
nificantly lower in the secondary peak, i.e., the spectrum in this 
phase bin is harder than in the other phase bins. The centroid en- 
ergy of the lOkeV feature does not vary with pulse phase while 
its depth weakens during the main peak. The equivalent width 
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Fig. 8. a: Difference of x 2 of the best fit of epoch I with frozen 
CRSF parameters (Fig. |2jc) and the x 2 of the best fit without 
CRSF. The centroid energy has been varied in steps of 1 keV. 
The red, dashed-dotted curve displays the y 1 of th e best fit when 
using tcrsf and ctcrsf from iHeindl et al.l ([2001). For the light 
blue, solid curve the parameters from this work have been used. 
The dark blue, dashed-doubledotted curve shows this difference 
when both, the width and depth (the upper limit has been set to 
3) of the CRSF are allowed to vary, i.e., the same situation as 
that of Fig. [2] The two additional minima in x 1 are discussed 
in Sect. 13.21 The horizontal dashed line indicates Ay 2 = 0, and 
the horizontal dotted line seperates the logarithmic and the lin- 
ear scaling in y-direction. The ve rtical dashed lines s hows the 
position of the CRSF derived in IHeindl et al.l ( 1200 lb and this 
work, b and c: results for the optical depth and the width of 
the CRSF, respectively, when allowing these parameters to vary. 
The plateau of the curve displayed in panel b at ~60-70keV is 
caused by the upper limit for tcrsf of 3. 



remains constant in the main and the secondary peak, while it 
shows a larger value in the off state. Finally, we note that there 
is neither evidence for changes of the Fe line flux nor for the 
equivalent width over the pulse phase. 

4. Results and conclusions 

4.1. Outburst series 

In this paper we analyzed quasi-simultaneous XRT, ISGRI, and 
PCA observations of two outbursts of XTE J1946+274 during a 
series in the second half of 2010. 

Before the onset of this outburst series, XTE J1946+274 was 
in a state of quiescence for almost one decade, i.e., about 20 or- 



7 



Sebastian Miiller et al.: The reawakening of the sleeping X-ray pulsar XTE J 1946+274 



bits of the neutron star. Even though the formation, structure, 
and dynamics of Be disks ar e even today far fr om being com- 
pletely understood (see, e.g., iDraper et alj|201 ll and references 
therein), the missing mass accretion onto the neutron star during 
that time is probably due to the absence of a sufficiently large Be 
disk during that time. A new outburst was only possible once the 
disk had been replenished and accretion could be triggered. 

This series shows similar behavior as the one observed in 
1998, e.g., two outbursts are observed per orbital period. The 
outbursts of XTE J 1946+274 could also be similar to tho se seen 
in GX 301-2, which have been extensively modeled bv iLeahvl 
d2002l) . These authors posit that an additionally stream of matter 
is flowing from the primary, and that a second outburst per orbit 
could be caused by the passage of the neutron star through this 
stream. These flux peaks of GX 301-2 occur near the apsides of 
this system. 

However, the outbursts of XTE J 1946 +274 do not clearly co- 
incide with the times of periastron and apastron passages of the 
neutron star (Fig.[T). This could be explained by 3 -dimensional 
simulations, which show that the disturbance of the Be disk by 
the gravitational field of the neutron star could lead to a strong 
asymmetric structure of the circumstellar material, which could 
also lead to multiple X-ray o utbursts during one orbital period 
(see, e.g jQkazaki et alJl201 ll) . 

Alternatively, the outbursts can be triggered by the neutron 
star passing through the Be disk due t o a misalignment o f the 
orbit and the Be star's equatorial plane ( Wi lson et al.ll2003h . 

Another possible explanation for the irregularity of the out- 
bursts is that they are triggered by density variations in the Be 
disk and not by orbital effects alone. This assumption could be 
verified by regular optical monitoring of the Balmer Ha line, 
which is an indicator fo r the presence of such a disk (see, e.g., 
iGrundstrom et al. 2007, and references therein). 

4.2. Spectroscopic results 

An absorbed Fermi-Dirac cutoff powerlaw together with an iron 
Ka fluorescence line, an iron line complex between 6 and 7 keV 
caused by Galactic ridge emission, a Gaussian like absorption 
feature around lOkeV, and a cyclotron line at ~25keV repro- 
duce the observations well in terms of xL*. We find time as 
well as pulse phase dependent variability of the continuum pa- 
rameters of XTE J1946+274. In the time resolved case, these 
changes might be caused by different accreting mechanisms de- 
pending on the mass transfer onto the neutron star. The periodi- 
cally changing line of sight with respect to the X-ray emitting re- 
gion, caused by the rotation of the neutron star, is likely respon- 
sible for the variabilities observed in the pulse phase resolved 
analysis. Furthermore, this rotation together with differences in 
the accretion geometries at the two magnetic poles of the neu- 
tron star might lead to the asymmetric pulse profile of the X-ray 
pulsar (see Fig.0. 

Studying the behavior of the spectral shape is problematic, 
because the spectral parameters photon index T, folding energy 
Zifoid, as well as the cutoff energy E cut , which describe the broad- 
band shape of the X-ray spectrum, show strong cross correla- 
tions to each other. To account for this, we calculated ratios of 
the spectra for both, the time resolved and the pulse phase re- 
solved studies. In the time resolved case, we find clear devi- 
ations from a constant ratio for epoch V. During epochs I-IV, 
the overall spectral shape seems to remain relatively constant. In 
the pulse phase resolved analysis, the spectrum of the secondary 
peak turns out to be significant harder than these from the main 
peak and the off state. Such properties of X-ray spectra are due 



to, e.g., the temperature of the visible part of the X-ray emitting 
region, in particular the accretion column. 

Another feature required to get a good fit is a source intrin- 
sic Gaussian iron Ka emission line at 6.4 keV. The disagree- 
ment in the strength of the required emission in this energy 
range between PCA and XRT in initial fits can be solved by in- 
ducing the emission of the galactic ridge. The strength of the 
Galactic ridge emission is about four times larger th an found 
in other regions of the Galaxy (lEbisawa et al.ll2007h . Such a 
difference is consistent, ho wever, with the typica l spatial varia- 
tions of the ridge em ission dYamauchi et al.ll2009t) . Furthermore, 
iKiihnel et al.l (1201 2h found an emission strength consistent with 
our result for GRO J1008-57. 

After taking into account the contribution due to Galactic 
ridge emission, the flux of the source intrinsic Fe Ka emission 
line is correlated with the X-ray flux Fx (s ee Fig. |7]>, a s also 
observed in other X-ray transients (see, e.g., IInoudll985l) . The 
Fe Ka line is significantly detected at source fluxes (7-15 keV) 
greater than about 0.035 photons s cirT 2 . At fluxes below this 
value, the source intrinsic line is consistent with zero. The equiv- 
alent width Wpe stays constant for high lumi nosities, i.e., epochs 
I— III. Reanalyzing the earlier RXTE data dHeindl et al.l 1200 lh 
with the continuum model employed here yields an equivalent 
width which is also in agreement with the result from the 1998 
outburst (W£ 9 * = 59 + 10 eV). The different values for W Fe for 
the fainter observations, where the Fe Ka line is only marginally 
detected, can be caused by the uncertainties of the Galactic ridge 
emission, which significantly contributes to the data at these en- 
ergies. In the pulse phase resolved analysis we find no evidence 
for variability of the line flux. This result indicates that it is emit- 
ted in a region large compared to the distance travelled by light 
during one pulse period. The equivalent width also shows no 
variation, which we would actually expect for a constant model 
component at varying flux levels. However, this can be explained 
by the quite large relative uncertainty of Wp e , which is on the 
same order of magnitude as the respective variations. 

The lOkeV feature is present during all epochs and pulse 
phases. It shows relatively constant results for line energy, width, 
and depth in the time resolved case. The reanalysis of the 1998 
data yields an equivalent width of W^ 9& = -250^™ eV, which 
also equals the results for the current outburst. While we find al- 
most no connection of the parameters of this feature with the lu- 
minosity, its optical depth and equivalent width varies with pulse 
phase. Wg is constant during the main and the secondary peak, 
but it strongly increases during the off state. This behavior indi- 
cates that the fractional amount of absorbed flux related to the 
continuum level is constant for the main and the secondary peak 
and changes in the off state. This behavior of the lOkeV feature 
could give rise to speculations about possible physical processes 
as, e.g., pulse phase dependent absorbing processes, producing 
such a feature. On the other hand, the residuals around lOkeV 
could also be due to the failing of the spectral broadband contin- 
uum models resulting in a wrong description of the data around 
these energies. More quantitative analyses of this feature, also 
from other sources where it occours in emission, are urgently 
needed to reveal the true nature of this enigmatic feature. 

4.3. Cyclotron resonance scattering feature 

We find w eak evidence for the presence of the CRSF, first dis- 
covered bv lHeindl et ail (120011) during the 1998 outburst series. 
The cyclotron line improved the model fit for the high signal to 
noise epoch I spectrum, where PCA and ISGRI data are avail- 
able. Later spectra are consistent with the presence of a CRSF 
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with unchanged parameters, however, the line is not formally 
detected in these observations because of their lower signal to 
noise ratio. We stress that due to its relative weakness, the inclu- 
sion of the CRSF in our fits does not affect our results for the 
continuum parameters or for the pulse phase resolved analysis. 

If the identification of the line is correct, then its depth in 
2010 was significantly lower than that measured in the 1998 
outburst (r^p = 0.09+° '° vs. t 1 ™ f = 0.33 + ° -jg. IHeindl et al.l 
l200l . Furthermore the centroid energy of the CRSF in 2010 
( £ c°rsf = 25.3+» 9 | ke V) is significantly lower than in 1998 

( £ crsf = 36 - 2+ o 7 keV. IHeindl et alj|200lh . Reanalyzing the ear- 
lier RXTE data as described above yiel ds CRSF parameter s 
which are consistent with those found by IHeindl et al.l (1200 ll) . 
Furthermore, the equivalent width in 1998 was also significantly 
larger than in the current data (H^rsf = (-2.1+^) x 10 3 eV). 

One reason why constraining the CRSF is so much eas- 
ier in the 1998 data is due to the fact that, while the expo- 
sure times of the 1998 and the 2010 data are comparable, the 
10-20keV X-ray flux of the 1998 observation is -40% higher 
than in the 2010 observation (~9.6 x 10~ 2 photons s _1 cm~ 2 vs. 
~6.9 x 10~ 2 photons s _1 cm 2 ). However, fits to the 2010 data in 
which we fix the CRSF parameters to their 1998 values do not re- 
sult in a satisfactory description of the data (x^. ed = 1.55). Could 
this indicate that the energy of the line varied? Assuming that 
the CRSF width and depth remained the same, we searched for 
a line similar to that found in the 1998 data. Varying the CRSF 
energy in steps of a 1 keV and refitting the continuum parame- 
ters did not yield any satisfactory results, with best-fit x^ ed > 1-5 
over the 20 to 60keV band (see Fig. [8] red curve). This result 
indicates that irrespective of the difference in signal to noise be- 
tween both data sets the CRSF must have varied between both 
outbursts, and a strong line as that seen in 1998 is not consistent 
with the data analyzed here. 

Taking the 2010 CRSF values at face value, a possible ex- 
planation for the difference between the 1998 and 2010 out- 
bursts could be the flux d ependence that is seen in s ome other 
CRSF sources (see, e.g.. ICaballero & Wilms! I20TTI for a re- 
cent discussion). XTE J1946+274 could therefore be a cyclotron 
source with an overall positive correlation betw een the X-ray 
flux a nd the CRSF energy, similar to Her X-l (IStaubert et al.l 
120071) . In models for the change in cyclotron line energy in 
neutron stars this change is generally interpreted as a change 
in height of the line producing region as the accretion column 
changes with mass accretion rate (and thus luminosity). Using 
the dipole approximation, and the assumption that in 1998 the 
CRSF was emitted from the neutron star's surface (r^s = 10 km 
and tons = 1.4 M Q ) yields a height difference of 1.2 km for the 
regions where the CRSF is generated. This result is in agree- 
ment with typical estimates of several kilometers for the heigh t 
of the accre ti on co lumns (see, e.g., [Basko & Sunvaevlll976l) . 
iBecker et al.l (1201 2h show that a positive correlation between 
flux and energy is possible in the luminosity range where the 
stopping in the accretion column is dominated by Coulomb 
braking. For cyclotron line energies around 30keV Coulomb 
braking is the dominating braking process in the luminosity 
range 1-5 x 10 37 ergs~'. For higher luminosities, radiation brak- 
ing dominates. XTE J1946+274's peak outburst luminosity of 
4.5 x 10 37 ergs~' is barely consistent with this range, however, 
given that the source seems to be located at the transition be- 
tween Coulomb and radiation braking and was brighter in 1998, 
the very large change in CRSF energy seems unlikely to be due 
to a pure mass accretion rate effect. 

In conclusion, while the possibility of a luminosity depen- 
dent CRSF is intriguing, the poor signal to noise ratio of the 



2010 data does not allow a definitive answer concerning the lu- 
minosity dependence of the CRSF. Further, longer, monitoring 
observations of XTE J 1946 +274 during its next outburst episode 
are urgently needed to resolve this question. 
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